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ABSTRACT 

We derive atmospheric parameters and lithium abundances for 671 stars and include our measure- 
ments in a literature compilation of 1381 dwarf and subgiant stars. First, a "lithium desert" in the 
effective temperature (T e g) versus lithium abundance (^4l0 plane is observed such that no stars with 
T c g ~ 6075 K and Au ~ 1.8 are found. We speculate that most of the stars on the low Au side of the 
desert have experienced a short-lived period of severe surface lithium destruction as main-sequence 
or subgiant stars. Next, we search for differences in the lithium content of thin-disk and thick-disk 
stars, but we find that internal processes have erased from the stellar photospheres their possibly 
different histories of lithium enrichment. Nevertheless, we note that the maximum lithium abundance 
of thick-disk stars is nearly constant from [Fe/H] = —1.0 to —0.1, at a value that is similar to that 
measured in very metal-poor halo stars (Ay\ — 2.2). Finally, differences in the lithium abundance 
distribution of known planet-host stars relative to otherwise ordinary stars appear when restricting 
the samples to narrow ranges of T e g or mass, but they are fully explained by age and metallicity 
biases. We confirm the lack of a connection between low lithium abundance and planets. However, 
we find that no low Au planet-hosts are found in the desert T e g window. Provided that subtle sample 
biases are not responsible for this observation, this suggests that the presence of gas giant planets 
inhibit the mechanism responsible for the lithium desert. 

Subject headings: stars: abundances — stars: evolution — stars: planetary systems 



1. INTRODUCTION 

The current lithium abundance in the solar photo- 
sphere is over two o rders of magnitude l ower than the 
meteoritic value (e.g., IAsplund et aLl feOOQ). implying sig- 
nificant lithium depletion since the solar system formed. 
For many decades, this observational fact has posed a se- 
vere problem for solar interior models. The temperatures 
needed for astration of lithium are found just below the 
base of the present-day solar convective zone. Therefore, 
according to standard models, there should be no surface 
lithium depletion during the main-sequence life of stars 
like the Sun (e.g.. iD'Antona fc Mazzitellilll9940 . Yet, in 
the solar neighborhood, lithium abundances have been 
found to vary by more than two orders of magnitude 
among stars with similar atmosphe ric parameters, in 
parti cular near the solar values (e.g.. lLambert fe Reddvl 
2004). Several mechanisms have been proposed to ex- 
plain the observed lithium abundances in solar-type 
stars, including rotational mixing (e.g.. | Pinsonneaultl 
120101 ) . gravity waves (e.g.. iCharbonnel fc TalonT l2005L 
and atomic diffusion (e.g., Mic haudl 11986( 1 . but a fully 
consistent picture remains elusive. 

Another important problem for models of lithium 
depletion is the observation of a lithium chasm in 
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the effective temperature (T c s) range between about 
6500 and 6850 K in stars i n young open clusters, par- 
ticularly the Hyades (e.g., iBoesgaard fc Tripiccol 119861 : 
iSoderblom et al.lfl993l ) . Within this narrow temperature 
range, whose center corresponds to about 1.4M in the 
case of the Hyades, stars have significantly lower lithium 
abundances than those immediately outside of that win- 
dow. That this so-called "lithium dip" is developed in 
main-sequence rather than pre-main sequence stars is 
demonstrated by the fact that it is not seen among stars 
of the youngest open clu sters (jBoesgaard et alJ " ll988t 
Balachan dran et al.ll201ll ). The lithium dip is, as ex- 
pected, seen in field main-sequence and s ubgia nt stars 
ilBalachandrarJ 119901 : lLambert fc Reddvl I2004D where 
the dip moves to l ower mass with lower metallicity 
(jBalachandranl [l995h . Attempts have been made also 
to explain the lithium dip in young open clus ters (e.g., 
iTalon fc Charbonnell H99cl iMendes et al.lll999fl . but the 
overall picture is still far from being understood com- 
pletely. 

The complex behavior of lithium abundances observed 
in FGK stars and the poorly constrained mechanisms of 
lithium depletion inside the stars make it difficult to use 
stellar lithium abundances to constrain Galacti c chem- 
ical evolution models (e.g.. iRomano et al"1ll999r i. Sites 
of lithium production have been identified, but their 
relative importance for the lithium enrichment of the 
interstellar medium is still somewhat debatable (e.g., 
iTravaglio etlll[200lh . To the best of our knowledge, 
there has not been a study of lithium abundance evolu- 
tion in different local disk stellar populations. 

Studies of samples of stars with lithium abundances 
homogeneously derived (within each study) have been 
published in the past few years. Analysis of these samples 
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has allowed a better understanding of lithium depletion, 
but an overall picture still needs to be drawn from the 
wealth of av ailable data. 

Recently, iRanrirez et al.1 (|2007t ) have extended their 
high resolution spectroscopic study of nearby stars from 
about 500 to 800 objects (Ramirez, Allende Prieto, & 
Lambert, in preparation, hereafter RAL). Their data are 
being used to study oxygen abundances in local stellar 
populations, but their spectra cover the 6708 A region, 
which contains the lithium doublet that is commonly 
used to derive lithium abundances. In this paper, we use 
RAL's dataset as a starting point to construct the largest 
available catalog of lithium abundances in nearby FGK 
dwarf and subgiant stars. The catalog is constructed 
using also data from the literature, but RAL's stellar pa- 
rameters and the lithium abundances derived here are 
used to normalize all data sets to a common scale. 

Given the nature of the RAL work, our sample is 
suited to study systematic differences in the chemical 
evolution of solar neighborhood stellar populations. At 
least two distinct populations of stars, with different 
spatial distributions and kinematical properties, co-exist 
in the Galactic disk, namely the th i n disk and th e 
thick disk (e.g., IGilmore fe Reidl [l98l ISoubiranl Il99l . 
Moreover, it is known that these two populations differ 
in chemical comp osition, particularly in the q-elements 
and oxygen (e. g.JBensbv et al.l[2005tlR,eddv et al.112001 
iRamirez et al.ll2007t ). Models of Galactic formation and 
evolution must be able to explain these chemical and 
kinematical differences. As mentioned before, we are not 
aware of any prior analysis of differences in lithium abun- 
dances between the thin and thick disks, which can be 
investigated by studying stars that have experienced the 
least amount of lithium depletion. RAL's data set in- 
cludes also spectra of a number of old halo stars, which 
we can use to shed new light on the problem of the pri- 
mordial, or cosmologica l, lithium abund a nce. 

Som e authors (e.g., 1 Gonzalez! 120081 : llsraelian et al.l 
2009) claim to have revealed a connection between stellar 
lithium abundance and the presence of exoplanets in Sun- 
like stars, although others have not found evidence fo r 
this correlation (e.g.. lRvanll2000t iBaumann et al.ll2010D . 
The mechanisms proposed to explain how the presence 
of planets could affect lithium depletion are planetary 
migration, which affects the evoluti on of the angu l ar mo - 
mcntum of the star, as sugg ested bv lCastro etaLl (|2008h . 
and interactions between the protoplanetary disk and the 
star, which determine the degree of differential rotation 
between the radiative core and the convective envelope, 
therefore having an imp ortant impact on rotational mix- 
ing fe.g.. lBou"vi cr 2008). In both cases, enhanced lithium 
depletion is expected for stars hosting planets. Lithium 
abundances in stellar photospheres could, therefore, be 
of interest also for exoplanet research. However, sam- 
ple biases must be clearly identified and taken into ac- 
count before drawing any strong conclusions from lithium 
abundance data sets. The large size of our sample allows 
us to control these biases and to investigate any possible 
underlying connections between the presence of exoplan- 
ets and stellar surface lithium abundance. 

Clearly, lithium abundances are very important for 
studies of stellar evolution and physics. By analyzing 
them, we hope to gain insights into not only stellar struc- 



ture and evolution, but also the chemical evolution of 
lithium in the interstellar medium and, arguably, the im- 
pact of planet formation on stellar interior evolution. 

2. SAMPLE, OBSERVATIONAL DATA, AND LITHIUM 
ABUNDANCE ANALYSIS 

The data set analyzed in this work consists of about 
700 high resolution {R = A/AA > 45,000), high signal- 
to-noise ratio (S/N > 100) spectra of 671 nearby (d < 
200 pc) FGK dwarf and subgiant stars taken from the 
RAL study. The high quality of these data allows a pre- 
cise determination of lithium abundances even in stars 
with severe lithium depletion. 

Most of the spectra were taken by us using the Tull 
coude spectrograph on the 2.7m Harlan J . Smith Tele- 
scope at McDonald Observatory (jTull et alJ li995). We 
also used FGK stars' spectra fro m the online UVES-VLT 
library by iBagnulo et al.l (|2003l ) and the library of spec- 
tra of nearby stars — the Spectrosco pic Survey of Stars in 
the So lar Neighborhood, S 4 N — by I Allende Prieto et al.l 
(2004), taken from both the Smith Telescope at McDon- 
ald Observatory and ESO's 1.5 m Telescope at La Silla. 
Some of RAL's spectra were not useful for lithium abun- 
dance determination due to relatively low S /N or incom- 
plete wavelength coverage. Because of the difficulty of 
combining spectral data from sources with different spec- 
tral resolutions, when we had multiple spectra for a star, 
instead of co-adding the spectra, we independently ana- 
lyzed them and averaged the abundances obtained. 

We determined the lithium abundance using spectral 
line-profile fitting of the region near the lithium reso- 
nance doublet at 6708 A, as in IBaumann et aTl 12010). 
We u sed the spectrum synthesis code MOOG ()Snedenl 
[1971 2010 versio r0), which assumes local thermody- 
namic equilibrium (LTE), while the stellar parameters 
used are from RAL and were derived as described below. 

The effective temperature (T c g) was obtained from 
the metallicity-depend ent color- T e g calibrations by 
iCasagrande et afl (|2010f ) . Most of our sample stars have 
ground-based measurements of (B — V) and (6 — y) col- 
ors avail able in the General Catal og of Photometric Data 
(G CPD. iMermilliod et al.l 11997ft . the Hipparcos Cata- 
log dPerrvman et al.lll997f). and the Geneva- Copenhage n 
Survey (jNordstrom et al.ll2004J : ICasagrande et al.ll2011l) . 
Stars fainter than about V = 6 also have reliable infrared 
(JHK) photometry available from the Two-Micron All- 
Sky Survey (2MASS ICutri et al.ll200l . A small fraction 
of our sample stars also have (RI)c photometry available 
in the GCPD. We used as many as possible of the color 
combinations that are available with these data to de- 
termine Tea, averaging different T e ff values derived from 
all available colors for each star (weighted by their er- 
rors, obtained by propagating the photometric uncertain- 
ties, the scatter of the color-T e ff calibrations, and con- 
sidering also the error in [Fe/H]). The surface gravity 
(log<?) was derived from mass estimates based on the- 
oretical isochrone analysis and trigonometric Hipparcos 
parall axes. The isochrone set used is that by iKim et al.1 
(2002f). Metallicities ([Fe/H]j3 were derived from equiv- 
alent width analysis of a large number of Fe I and Fe II 

6 http:/ /www. as. utexas.edu/~chris/moog. html 

7 In this work we use the standard definitions: [X/Y] = 
log(JVx/iVV) - log(Nx/N Y )e, and A x = log(Nx/Nn) + 12, where 
A^x is the number density of element X. 
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lines. Microturbulente velocities (v t ) were determined in 
the standard manner, i.e., by looking at the correlation 
between iron abundance and line strength of Fe I lines in 
each spectrum, and iteratively increasing or decreasing 
the v t value until the slope of that correlation is close to 
zero. 

The iron abundances inferred separately from the neu- 
tral and singly ionized species of iron agree satisfacto- 
rily. We find that their mean difference (Fell-Fei) is 
0.01 ± 0.06 dex, implying that ionization balance is sat- 
isfied for the majority of the sample stars. The aver- 
age internal errors in T e ff, \ogg, and [Fe/H] are 53 K, 

0. 08 dex, and 0.05 dex, respectively. The mean value 
and error in [Fe/H] were obtained from the average and 
line-by-line scatter of Fe I and Fe n lines combined, re- 
spectively. Using all iron lines together for these calcu- 
lations is justified by the fact that ionization balance is 
satisfied. The procedure described above was iterated 
several times because the stellar parameters are interde- 
pendent. We made sure that convergence to a final result 
was achieved for each star. 

The model stellar atmospheres employed are from 
the MARCS grid with standard chemical composition, 

1. e., taking into account a-element enhancement at low 
[Fe/H]0 Gaussian macroturb ulent velocities wer e deter - 
mined using the formulae by iValenti fe Fischer! (|2005l) . 
defined as the lower envelope of the macroturbulence as 
a function of T c r derived by ignoring rotation. Rota- 
tional broadening was determined using line profile fit- 
ting of nearby Fei lines, after applying instrumental and 
macroturbulent broadening. 

Errors in lithium abundances due to uncertainties in 
stellar parameters were determined as follows. For each 
star, we generated a synthetic spectrum of the 6708 A 
region using our derived stellar parameters and lithium 
abundance, without external broadening convolution. To 
calculate the error in lithium abundance due to the un- 
certainty in T e fi, we increased the temperature of the 
model to T G ff + AT G g- where AT gf is the error in effec- 
tive temperature according to RAL, and computed the 
lithium abundance necessary to match the first profile 
(that with our derived lithium abundance) . We repeated 
the exercise for T e g — AT c g. The average of the dif- 
ferences between these new lithium abundances and the 
original one was adopted as the error due to the T e ff 
uncertainty. A similar procedure was used to measure 
the smaller errors due to \ogg and [Fe/H] uncertainties. 
The three errors were then added in quadrature. This 
is a reasonable approximation considering that the at- 
mospheric parameters were derived by RAL using nearly 
independent techniques 

For the majority of our sample stars, the error in 
lithium abundance due to the synthesis fitting proce- 
dure itself (i.e., the error due to observational noise) was 
found to be negligible compared to the error due to stel- 
lar atmosphere uncertainties. A few noisy spectra were 

8 http://marcs.astro.uu.se 

9 The stellar parameters derived by RAL are not fully indepen- 
dent, but their correlation is weak and their impact on the errors is 
smaller than the internal uncertainties. For example, for a Sun-like 
star, a change of 50 K in T e g corresponds to 0.03 dex for [Fe/H] (if 
inferred from Fei lines, but smaller if derived from Fell). Propa- 
gating this [Fe/H] shift to re-determine T c ff leads to a change of 
8K, which is a small fraction of the original T c g error. 



Table 1 

Stellar Parameters and Lithium Abundances Determined in 
This Work. 1 



HIP 



T eff (K) log g [cgs] 



[Fe/H] 



171 
348 
394 
475 
493 
522 
530 
544 
656 
768 



5510 ± 66 
5746 ± 55 
5635 ± 37 
5836 ± 72 
5960 ± 44 
6251 ± 44 
5866 ± 40 
5458 ± 40 
5805 ± 39 
5681 ± 31 



4.46 ±0.01 
4.38 ±0.06 
3.78 ± 0.07 
4.35 ± 0.05 
4.41 ± 0.03 
4.21 ±0.02 
3.90 ± 0.05 
4.52 ±0.02 
3.82 ±0.06 
3.99 ±0.04 



-0.76 ±0.06 
-0.19 ±0.04 
-0.48 ±0.04 
-0.34 ±0.06 
-0.20 ±0.04 

0.05 ± 0.05 
-0.01 ± 0.05 

0.14 ±0.06 
-0.24 ±0.06 
-0.01 ± 0.06 



< 0.35 
0.82 ±0.04 
1.98 ±0.04 
1.63 ±0.05 
2.32 ±0.04 
2.84 ± 0.04 
0.70 ± 0.03 
2.39 ±0.04 
2.53 ±0.04 
1.46 ±0.02 



1 Table [T] is published in its entirety in the electronic edition of 
the Astrophysical Journal. A portion is shown here for guid- 
ance regarding its form and content. 

the exception. They were assigned a conservative error 
of 0.15 dex. This is larger than any error determined 
for our sample stars for which lithium abundances were 
measured. In cases where the lithium doublet was not 
detected, we were still able to calculate an upper limit 
to the lithium abundance based on the stellar parame- 
ters of the star and the local signal-to-noise ratio of the 
spectra (these values correspond to 2 a detection limits) . 
Our lithium abundance results are given in Table HE] 



3. A CATALOG OF LITHIUM ABUNDANCES 
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We supplemented o ur measurements 
lithium abund anc es fro m lLambert fe Reddvl ( 2004 ) 
i Luck fe Heiterl (120061): llsraelian et all ( 2009 ) 
Baumann et al l (120101) : iGhezzi et"aTf (201C) 
Gonzalez et all (120101): and iTakeda et all (|2010T) . 
lLambert fe Reddvl (120041) inc luded a com pilatio n of 
data fromlBalachandraTjl990l ): IChen et all 1)20011 ): and 
iReddy et all ( 2003[ ). which they ensured to be consistent 
with regards to the stellar parameters and lithium 
abundances. All of these works used high quality 
data and their results are, internally, analyzed in a 
homogeneous fashion. These previously published data 
were normalized to ours by computing and subtracting 
the mean offsets in T c g, \ogg, [Fe/H], and lithium 
abundance, to compensate for systematic differences in 
stellar parameters which may affect our mass and age 
determinations. Stars in common between this work and 
the previously published studies were used to calculate 
these offsets, which are given in Table [21 

Instead of deriving lithium abundance offsets based 
on the mean differences in stellar parameters, we used 
the mean offsets in lithium abundances themselves to 
put all lithium abundance measurements into the same 
scale. This is because the lithium abundances depend not 
only on the stellar parameters adopted by each group, 
but also on the atomic data employed in the spectral 
line synthesis of the 6708 A lithium doublet. Upper 
limits to lithium abundances were excluded from the 

10 RAL's data set also includes 9 solar spectra, either from 
(reflected sunlight) asteroid or skylight observations. The aver- 
age LTE lithium abundance for these spectra is 1.02 dex, with a 
standard error of 0.02 dex. This value is in good agreement with 
other recent determ inations of the solar lithium abundance (e.g., 
lAsplund etani2009h . 
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Table 2 

Offsets in Stellar Parameters and Lithium Abundances Relative to This Work. 1 



Samplc a 


AT cff (K) 


Alogg 


A[Fe/H] 


AA Li 


JV b 


N to t c 


BIO 


40 ±48 


0.05 ± 0.06 


0.03 ±0.03 


0.07 ± 0.09 


69 


117 


G10 


-4 ±68 


-0.02 ± 0.04 


0.04 ±0.04 


-0.13 ±0.12 


44 


152 


109 


26 ±32 


0.05 ±0.05 


0.02 ± 0.02 


0.06 ±0.04 


11 


80 


LH06 


50 ±85 


0.04 ±0.19 


0.01 ±0.06 


0.01 ±0.11 


76 


194 


T10 


-12 ± 54 


-0.03 ±0.09 


0.03 ± 0.04 


-0.03 ±0.10 


26 


117 


GhlO 


30 ± 71 


-0.06 ±0.14 


-0.00 ±0.06 


0.01 ±0.08 


38 


259 


LR04 


-83 ± 63 


-0.01 ±0.11 


-0.02 ±0.05 


-0.07 ±0.08 


212 


449 



1 The errors quoted for the A values correspond to the 1 a standard deviation of the 

differences. 

a Sample: B10 = [Baumann etall J2010l), G 1 = fGonzalez efall <201Ch, 109 = 
Israelian ct al. (2009), LH06 = Luck & Hciter (2006), T10 = Takcda ct al. (2010), 
GhlO = IGhezzi et all | |2010|) . LR04 = ILambert fc Reddvl pOOlfr 



b N is the number of stars in common with this work. 
c JVtot is the number of stars in our final catalog taken from each study. 



co mputation of the off sets. ILambert fc Reddvl (j2004f ) 
and lTakeda et al. (2010) applied non-LTE corrections to 
their lithium abundances, which we undid in order to pre- 
vent different non-LTE prescriptions from systematically 
biasing our data, i.e., we recovered their LTE lithium 
abundances. 

When a star was available in more than one source, we 
calculated a weighted average for each of the parameters 
and lithium abundance. We adopted as weights the in- 
verse vari ances given in each o f the sources. The errors 
quoted by Israelian et all (|2009f > are very small compared 
to those from other sources. Such small errors are prob- 
ably reasonable for internal analysis of their sample, but 
systematic errors will dominate when comparing to other 
works. In orde r to ay oid giving too much weight to the 
llsraelian et al.l (|2009ft data (no source should be given 
very high preference in this type of compilation), we in- 
creased their error bars by a factor of 3 in T e ff and log g 
and a factor of 4 in [Fc/H] (retaining the relative size of 
the errors equal, therefore giving more weight to stars 
which were analyzed using higher quality data, or stars 
that are more similar to the Sun, where systematic errors 
are smaller). These multiplicative factors were estimated 
from the comparison of our data to theirs. The star-to- 
star scatter in the differences in the parameters for stars 
in common could be made compatible with the quoted 
errors if increased by tho s e fact ors. 

F or ILambert fc Reddvl ()2004D . iLuck fc Heiterl (2006), 
and lTakeda et all (|2oTojT we adopted errors of 50 K, 0.07 
(log g), and 0.05 ([Fc/H]) for all stars. Although more re- 
alistic errors, generally including systematics, of about 
100 K, 0.1 dex, and 0.1 dex, respectively, are typically 
quoted in the literature, we adopted these error bars to 
give a reasonable weight to these stellar parameter and 
lithium abundance measurements when computing aver- 
age values within our compilation. 

Our procedure of increasing errors as described in the 
previous paragraphs may appear somewhat arbitrary, al- 
though we have given a quantitative justification for it. 
Nevertheless, we note that not weighting the averages 
gives very similar results. Indeed, for the 475 stars which 
have data from more than one source, the difference in 
T e ff between the weighted and non-weighted values is 
± 12 K, that for log 5 is 0.1 ± 0.2, and for [Fe/H] we 
find —0.001 ± 0.011. These small mean differences have 
only a marginal effect on our subsequent analysis. 



In our compilation, when lithium upper limits were 
available from more than one source, we adopted the 
lowest upper limit found, which typically corresponds to 
the case where the highest quality data were used. 

To correct for non-LTE effects on the 6708 A lithium 
doublet formation and abundance derivatio n, we used 
the no n-LTE tables and interpolation code bv lLind et all 
( 2009). For a number of stars, the parameters (especially 
microturbulent velocity) were outside of the range cov- 
ered by these tables. For these stars, we used the nearest 
node in the stellar parameter grid of the non-LTE tables 
to determine the size of the correction. These extrapo- 
lations are not risky because of the asymptotic nature of 
the non-LTE correction s in these extreme s of the stellar 
parameter grid used by iLind et al.l (|2009f ) . For most of 
our sample stars, the non-LTE corrections are of order 
+0.02 dex, and more important (~ —0.05 dex) only for 
stars with very high lithium abundance (A-^i 2.5). 

Stellar masses and ages wer e derived as in RAL. 
The method is also de scrib ed in | Baumann et all (|2010[) ; 
iMelendez et~ai1 (|2012l ); and lChaname fc Ramfrel (|2012li . 
In summary, each star is placed on the HR diagram us- 
ing its T c ff and absolute magnitude. The latter was com- 
puted usi ng parallaxes from t he new reduction of Hippar- 
cos data (jvan Leeuwenll2007r ). All theor etical isochrone 
points from the Yonsei-Yale grid (e.g., lYi et all 120011 : 
iKim e t al. 2002) within a radius of 3 a from the observed 
parameters were then used to compute mass and age 
probability distribution functions, from which the most 
probable parameters and 1 u-like errors were computed. 
The same procedure was used to determine trigonomet- 
ric surface gravities. These logg values are internally 
consistent with the T e g and [Fe/H] parameters adopted 
for each star in this work, and we prefer them to those 
derived as the average of published values. 

Isochrone ages are unreliable for young stars. These 
objects, however, tend to be active and their rotation 
periods can be measured from the modulation due to 
starspots on photometric light curves and/or chromo- 
spheric emission data. Rotation periods can then be 
used to estimate stellar ages more precisely than using 
isochrones, particularly for young stars, because stars 
with convective e nvelopes slow d own in a predictable way 
as they age (e.g., [Barnes 20071). Therefore, if available, 
we used measured rotation periods (P r ot) and the age- 
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Table 3 

Catalog of Stellar Parameters and Lithium Abundances. 



HIP 


T eff (K) 


log 3 


[Fe/H] 




Mass Mq 


Age (Gyr) 


Pop. a 


Planets 


Source b 


171 


5449 ± 77 


4.46 ± 0.01 


-0.84 ±0.08 


< -0.51 


0.76+°°} 

— U.U1 


6.3 ± 1.4 


rot] 






TW+LH06 


348 


5740 ± 50 


4.38 ± 0.06 


-0.17 ±0.04 


0.85 ±0.07 


Q4+0 03 
u - 3 -0.02 


7 6fi+ 2,35 
' - DD -2.33 


[iso] 


thin 


no 


TW+B10 


394 


5636 ± 50 


3.76 ± 0.07 


-0.49 ± 0.04 


1.98 ±0.02 


1 10+ 015 
1 -- LU -o.os 


5.00i? ; 2 « 


[iso] 




no 


TW+LR04 


475 


5836 ± 72 


4.36 ± 0.05 


-0.34 ±0.06 


1.63 ±0.05 


92+ 03 
u ' 3 -0.03 




[iso] 


thin 


no 


TW 


493 


5943 ± 50 


4.38 ±0.03 


-0.23 ±0.04 


2.31 ±0.02 


n QQ+ 03 
u - yM -0.02 


°- ' ' -2.00 


[iso] 


thin 




TW+LR04 


522 


6266 ± 50 


4.21 ± 0.02 


0.04 ±0.04 


2.84 ±0.02 


1 24+ - 01 
1,z -0.01 


2 90+ ' 42 
z - 9u -0 .24 


[iso] 


thin 


yes 


TW+GhlO 


530 


5868 ± 50 


3.89 ± 0.05 


-0.05 ±0.04 


0.70 ±0.03 


1 20+ - 08 
1,zu -0.03 


4 on+ - 75 
^• 9U -0.4S 


[iso] 


thin 


no 


TW+LR04 


544 


5465 ± 50 


4.53 ±0.03 


0.14 ±0.04 


2.38 ±0.04 


96 +a02 
u - yo -0.03 


0.2 ±0.1 


rot] 


thin 


no 


TW+LH06 


656 


5819 ± 50 


3.82 ± 0.06 


-0.25 ±0.04 


2.53 ±0.02 


-, 9 c+0.09 


3 S9+1-53 
°- oy -0.33 


[iso] 






TW+LR04 


699 


6179 ± 63 


4.13 ±0.04 


-0.30 ±0.05 


2.33 ±0.10 


i ii+0.04 

± - i± -o.oi 


5 50+ ' 29 
°- ou -1.07 


[iso] 


thin 


no 


LR04 



1 Table [3] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

a Population: thin-disk, thick-disk, or halo, according to the star's kinematics. If not available, the kinematic criterion adopted 
in this work implies an ambiguous or unreliable association. 

b Source: TW = This W ork, B1 =IBaumann et al.l <201Ch . G 10 =IGonzalez et al.l puTOT) , 10 9 =IIsraelian et al. | ([2009^1 . LH06 
= ILuck fc Heiten (2006h . T10 = ITakeda et al.l H20101) . GhlO = IGhezzi et al.l l|2010n . LR04 = ILambert fc Reddvl pfjM l. 



P ro t relation by iBarned (|2007T > to estimate the ages of 
stars. Measure ments of rotat i on per i od adopted in this 
work are f rom iGaidos et all (120001) : iStrassmeier et al.l 
(pOOOh : and lPizzolato et alJ (|2003l l 

All the stars in our sample have kinematical data from 
RAL, namely radial velocities and Galactic space veloc- 
ities inferred from them and the stars' parallaxes and 
proper motions. From these data we compute a kine- 
matical probability that a star belongs to the thin disk 
(Pi), thick disk (P2), and halo (P3). The formulae used 
to calculate thes e proba biliti es are given in, f or exa mple, 
iMishenina et alJ ()2004D and iRaimrez et~aLl (|2007t ). To 
assign a given star to one of these grou ps, we adopt a 
criteri on stronger than that employed by Ramirez et al. 
(2007), which is Pi > 0.7 for thin-disk stars and P 2 > 0.7 
for thick-disk stars. In addition to those conditions, here 
we require P1/P2 > 10 for thin-disk stars a nd P2/P1 > 10 
for th ick-disk members, as suggested by iBensbv et al.1 
(2005). For halo stars we adopt P 3 > 0.5. With these 
constraints we find in our sample 898 thin-disk stars, 144 
thick-disk stars, and 43 halo members. The rest of our 
targets (296 stars) have ambiguous or uncertain popula- 
tion membership. 

The stars in our catalog are also flagged as known 
planet-hosts (Planets= "yes" ) , non planet- hosts (Plan- 
ets="no"), and stars whose planetary status is un- 
known (Planets=" ••■ "). The non planet-host 
stars have had their radial velocities monitored but 
no planets have been detected ar o und t he m yet. We 
used ta bles by |Fischer fc Valentil (120051): iSousa et al.1 
(120081): lIsTaelian et al.1 (120091): IGhezzi et al.1 J 20M: 
iBrueamver et al.l (|2011l) . and lWright et al.l(|2011[ ) to sort 
the stars into these categories. We found that 165 stars 
in our catalog are known to have planets while 360 stars 
are, as of this writing, deemed to have no sub-stellar mass 
companions. 

Our catalog of stellar parameters, lithium abundances, 
and derived masses and ages is given in Table [3] The 
errors for the T c ff and [Fe/H] values given in Table |3] 
correspond to the sample variance in cases where data 
from more than one source were available. However, 



we adopted minimum errors of 50 K in T c g and 0.04 dex 
in [Fe/H] to prevent very small unrealistic errors aris- 
ing from coincidental agreement of the statistically few 
sources available. The mean errors in T e ff, logg, and 
[Fe/H] given in this Table are 59 K, 0.08 dex, and 
0.05 dex, respectively, while that for the lithium abun- 
dances (excluding the upper limits) is 0.08 dex. In Ta- 
ble [3] we also provide population and planet flags, as de- 
scribed in the previous paragraphs, and the source(s) 
from which the data were taken. 

4. DISCUSSION 

4.1. The Lithium Desert 

Figure Q] shows our lithium abundances and effective 
temperatures. Overall, lithium is less abundant in cooler, 
less massive stars with larger convective envelopes. The 
two stars with T e g < 5500 K and high lithium abundance 
(A u > 2.8) are HIP 13402 and HIP 30034. Both stars are 
known to be very young, which means that their surface 
lithium abundances have not yet been depleted by inter- 
nal p rocesses. HIP 13402 has an age of a few million years 
(e.g., ICavrel de Strobel fc Cavrell|1989l) while HIP 3 0034 



is a pre-main sequence star (e.g.. iScholz et aLll2007t ). 

One of the most striking features o f our lithium abun- 
dance compilation, first suggested bv lChen et al.1 (|2001|h 
is seen clearly in Figure [TJ Between T D s ~ 5950 K and 
6100 K, the stars appear to separate into two groups 
of high and low lithium abundance, creating a "lithium 
desert," a region of the T e s-Au plane that seems to not 
contain any stars even though it is surrounded by them. 
This region corresponds to the area inside the polygon 
plotted in Figure [T] The lithium desert likely extends to 
slightly cooler and warmer effective temperatures, but 
its boundaries are blurred by our T e ff uncertainties. The 
extent of the lithium desert in both T c g and A^, and the 
abundance of stars immediately outside of its boundaries 
make it highly unlikely that this feature is a product of 
our sample selection and/or observational errors. 

Since many, but not all, of the stars in the lower part 
of the lithium desert have only upper limits for their 
lithium abundances, i.e., their A-^i values are based on 
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Figure 1. Lithium abundance as a function of effective tempera- 
ture for all stars in our catalog. Downward arrows represent upper 
limits to the lithium abundance. The polygon bounds the ap- 
proximate extent of the lithium desert and it is drawn here and in 
subsequent figures only to guide the eye. An error bar representing 
the average internal errors is shown at the bottom right. 

the non-detection of the 6708 A lithium doublet, the pos- 
itive slope of the Au versus T e g relation for those stars 
simply reflects the detection limit of the lithium feature 
for spectra with signal-to-noise ratio typical of that of 
our sample stars and those taken from the literature. 
The maximum value possible for a given star with 
no lithium doublet detection depends on T c g because 
the minimum line equivalent width detectable, which is 
a constant assuming all spectra have the same signal- 
to-noise ratio, requires a higher lithium abundance for 
higher T e g values. This probably explains the shape of 
the lower envelope of the A^ versus T e & relation plot- 
ted in Figure [TJ not only around the lithium desert, but 
also over the entire T e g range covered by the stars in our 
compilation. 

In order to investigate the nature of the lithium desert, 
as well as other patterns of lithium depletion across the 
stellar parameter space covered by our sample, in Fig- 
ure [5] we show versus T e g trends for stars of a given 
mass. Moreover, within each mass group, we divide the 
stars into two groups of metallicity, a "metal-rich" group 
of [Fe/H] > —0.1 (filled circles) and a "metal-poor" 
group of [Fe/H] < —0.1 (open circles). 

As expected, there is a strong correlation between T e g 
and mass - all our sample stars are on the main sequence 
or are slightly evolved. The well-known effect that at any 
given mass lower metallicity stars have higher surface 
temperatures is also obvious in Figure [2] 

Stars with masses between 0.70 and 1.35 M© can be 
found on the high lithium abundance side of the lithium 
desert, but, with very few exceptions, stars on the low 
lithium abundance side of the desert appear to have M > 
1.10 Mq. Moreover, at any given mass, stars in the lower 
part of the desert tend to be more metal poor than those 
with higher lithium abundances, an effect that is clearly 
seen in the case of 1.30 Mq stars. Due to their warmer 
effective temperatures, most stars with M > 1.4 M© are 
found beyond the lithium desert. 

Lithium abundance versus age trends for stars of a 



given mass (and metallicity) are shown in Figure OS Only 
stars with reliable ages (age/error > 3) are plotted there. 
For stars more massive than M — 0.8 M© there is some 
evidence of main-sequence lithium depletion, but the 
time-scales seem to depend not only on mass but also 
on the stellar metallicity. For M < 1.1 Mq , metal- 
rich stars seem to deplete lithium quicker than metal- 
poor stars, a feature that is particularly dramatic in the 
1.00 ± 0.05 Mq group. Interestingly, this pattern ap- 
pears to invert at higher masses, although the trends 
become less clearly defined. This picture is further com- 
plicated by the fact that the metallicity dependence of 
main-sequence lithium depletion seems to disappear for 
stars with M > 1.35 Mq. 

Figures [2] and [3] show that, in addition to the stellar 
mass, metallicity is a key parameter to understand sur- 
face lithium depletion. This is not entirely surprising, 
since the properties of the convective zone, in particu- 
lar the temperature at the b ase, have a high sensitiv- 
ity to both parameters (e.g., iPinsonneault et ail 120011 : 
iTheado fc Vauclairl [20121) . 

Although the morphology of the lithium desert is ac- 
centuated by the combination of low and high mass stars 
in Figure [U it can be still detected in plots showing stars 
within a narrow mass range (Figure [5]). From M = 1.10 
to 1.30 Mq, stars within a mass range of 0.05 M© of a 
given mass are found both below and above the desert. 
Figure [3] does not show a sudden drop in the surface 
lithium abundance at a given age for stars in that mass 
range, although errors in our age determination might be 
preventing us from finding such trend. In these groups, 
however, and in particular in the 1.10 Mq group, stars 
appear to separate into high and low lithium abundance 
subgroups in a discontinuous manner. It is also clear 
that there are no young stars (age < 2 Gyr) in the lower 
part of the lithium desert while those that appear rel- 
atively young (age ~ 3 Gyr) are also the most massive 
ones. Since the main-sequence lifetime of a 1.1 M© star 
is about 7.5 Gyr, the mechanism responsible for the big 
drop in surface lithium abundance, and hence the lithium 
desert, must act while the star is on the main sequence or 
as a subgiant star, and most likely not in the first 2 Gyr. 
We cannot rule out the possibility that this process is 
very short lived, as the dramatic, complete lack of stars 
i n the desert s uggest s. 

iChen et al.l (|2001[ ) argued that the stars on the low 
lithium abundance side of the desert are evolved lithium 
dip stars. They found a very strong and tight correlation 
between mass and metallicity for these low lithium abun- 
dance stars, which they used to support their idea, be- 
cause the mass of the lithium dip stars see ms to correlate 
with the metallicity of the cluster (e.g., Balachandran 
1995). However, using their data (i.e., their masses and 
metallicities) , we find a similarly tight correlation for 
their high lithium abundance stars, implying that the 
mass- [Fe/H] relation is inherent to the data set and sam- 
ple selection, rather than a property of groups of stars 
with different lithium abundance patterns. 

The mass-[Fe/H] correlation described above could also 
be due to systematic biases in the determination of 
masses. Using our data, we also find a correlation be- 
tween mass and [Fe/H] for the stars on the low lithium 
abundance side of t he desert, but the d ispersion is larger 
than that found by I Chen et al.1 (|2001l ). In any case, we 
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Figure 2. Lithium abundance as a function of effective temperature for stars of a given mass and metallicity. In each panel, stars have a 
mass within 0.05 Mq from the number shown on the top left. Filled (open) circles correspond to stars more metal-rich (metal-poor) than 
[Fe/H] = —0.1. Downward arrows represent upper limits to the lithium abundance. The polygon shows the location of the lithium desert. 



find a similar behavior for the stars with high lithium 
abundance. Another possible reason for this correla- 
tion could be simply that the T e g range of the desert 
is so narrow that the metallicity effect on the location 
of the main-sequence, which becomes redder as [Fe/H] 
increases, essentially projects into a mass dependence. 
Thus, our analysis suggests that the lithium desert can- 
not be explained only with evolved lithium dip stars. 
A possibly short-lived mechanism operating during the 
main sequence or subgiant phase of stellar evolution is 
likely responsible for this feature. 

4.2. Galactic Chemical Evolution of Lithium 

4.2.1. The Lithium Problem 

A conventional procedure for probing the Galactic 
chemical evolution (GCE) of elements is to establish 
the run of [X/H] or [X/Fe] with [Fe/H] from spectro- 
scopic analyses of samples of stars such as those in the 
solar neighborhood. In the case of lithium, since sur- 
face depletions are common and may be severe, the cus- 



tom has been to plot Au versus [Fe/H] and to adopt 
the upper envelope of the distribution as defining the 
growth of lithium abundance with the iron abundance. 
Then, this upper envelope is confronted with theories 
of lit hium nucleosynthesis and evolution of the Galaxy 
(e.g IRvan et al.lteOOltrTravariio et alJ200H:lAiibes et all 
12002 ). 

Our catalog of lithium abundances, if analyzed in the 
above fashion, would confirm previous attempts to map 
the GCE of lithium. Often overlooked in these attempted 
mappings is the question of lithium depletion in those 
stars defining the upper envelope. This issue is most 
obviously a concern at low metallicities where the up- 
per envelope i s defined by the so-called Spite plateau 
(jSpite & Spite] I1982D . for which A Li ~ 2.2. Now that 
observations from the Wilkinson Microwave Anisotropy 
Probe ( WMAP) have unequivocally defined the baryon- 
to-photon ratio, standard Big Bang nucleosynthesis pre- 
dicts a primordial li t hium abundanc e of Au ~ 2.7 (e.g., 
iRichard et all 120051 : iSteigmanl [2010h . and this has led 
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Figure 3. Lithium abundance as a function of age for stars of a given mass and metallicity. In each panel, stars have a mass within 
0.05 Mq from the number shown on the top right. Filled (open) circles correspond to stars more metal-rich (metal-poor) than [Fe/H] = —0.1. 
Downward arrows represent upper limits to the lithium abundance. 



to the so-called lithium problem (e.g., Fields! [201 lh . for 
which two broad solutions are in contention: 1) have 
the stars on the Spite plateau depleted their atmospheric 
lithium abundance by about 0.5 dex? or 2) do the stan- 
dard Big Bang nucleosynthesis calculations need recon- 
sideration on account of revisions to the standard mod- 
els of particle physics and/or cosmology? Given that 
the present (meteoritic) lithium abundance is Au ~ 3.3, 
an estimate free from uncertainty about stellar deple- 
tion effects, a change in the pre-Galactic abundance from 
— 2.2 to 2.7 profoundly reduces the inferred growth 
rate of lithium in the Galaxy. 

Depletion in metal-poor low mass stars would imply 
depletion may also occur in the less metal-poor and some- 
what higher mass stars which define the start of the up- 
per envelope's transition from the Spite plateau to higher 
lithium abundances in disk stars. Mass estimates for 
stars in our catalog show that the maximum mass of 
stars populating the Au versus [Fe/H] upp er envelope 
decreases with decreasing [Fe/H] (see also IChen et all 



12003 and ITambert k Reddvl[200l . from about 1.6 M 
at [Fe/H] = to 1.0 M at [Fe/H] = -1. At the lower 
masses, for example M < 1.3M Q , lithium is depleted in 
main-sequence stars with depletion increasing with de- 
creasing mass and increasing age, as clearly shown by 
the lithium a bundance reported f or ope n clusters — see, 
for example, Sestito k Randich (2005]), who, however, 
argue that depletion becomes ineffective for stars older 
than about 1-2 Gyr. Thus, the A Li versus [Fe/H] upper 
envelope most likely does not represent the growth of 
Galactic lithium abundance with [Fe/H]. The observed 
lithium abundances should be revised upwards by an 
amount that increases with decreasing [Fe/H]. 

4.2.2. Lithium in the Thin and Thick Disks 

It is well known that, at least in the solar neighbor- 
hood, thin- and thick -disk stars have different chemi- 
cal c omposition (e.g., iBensbv et all 2005; R eddv et al.l 
2006]). Thus, one expects that thin and thick disks may 
have different patterns for their lithium versus iron abun- 
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Figure 4. Lithium abundance as a function of [Fe/H] for thin-disk 
(upper panel) and thick-disk (lower panel) stars. Lines are drawn 
"by hand" to make the upper envelope differences clear. Solid 
(dashed) lines trace the upper envelope of the thin-disk (thick- 
disk) star data. Upper limits to lithium abundances based on the 
non detection of the lithium 6708 A doublet are not shown. 

dance evolution. Yet, to the best of our knowledge, no 
attempts have been made to find these different patterns. 
Our sample is of sufficient size and diversity that an at- 
tempt may now be made. Unfortunately, unlike the ele- 
ments which have been used to define the thin/thick disk 
differences in composition, lithium is subject, as noted 
above, to internal depletion; different Au versus [Fe/H] 
patterns for thin and thick disk may, therefore, reflect 
depletion dependencies rather than different histories of 
lithium nucleosynthesis. 

Lithium abundance against metallicity is plotted in 
Figure S] separately for thin- and thick-disk starsQ Our 
data clearly indicate different trends for the two disks. 
The full sample of thin-disk stars (Figure [H upper panel) 
shows an apparent enrichment with increasing metallic- 
ity, contrary to the case of thick-disk stars (Figure [4] 
lower panel) which appear to have a nearly constant max- 
imum lithium abundance. Interestingly, this maximum 
lithium abundance is similar to that of the Spite plateau 
exhibited by very metal-poor stars. 

However, the thin/thick disk difference shown in Fig- 

11 The thick-disk star with the highest lithium abundance in the 
lower panel of Figure Q] is HIP 10856 (A Li = 2.85 ± 0.04, [Fe/H] = 
+0.18 ± 0.06). This object has a high radial velocity and a high 
proper motion, which make its Galactic space velocities typical of a 
thick-disk star. Its atmospheric parameters, however, suggest that 
this object is younger than about 4Gyr. We note also that this 
star has a level of chromospheric emission above that of most stars 
of similar color (Isaacson & Fischer 2010), confirming its relatively 
young age. Moreover, its derived mass (1.18 ± 0.03 Mq) implies a 
lifetime of about 6 Gyr. Thus, instead of being a true thick-disk 
member, it seems more likely that the orbit of this star has been 
excited by a collision. If that were true, the highest Ali point in 
the bottom panel of Figure [4] should be in the top panel instead, 
where it would fit the overall trend perfectly. 
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Figure 5. Lithium abundance as a function of [Fe/H] for thin- 
disk (upper three panels) and thick-disk (lower panel) stars of mass 
M < 1.1 Mq. In the upper three panels, only thin-disk stars older 
than 2, 5, and 8 Gyr, respectively, are shown. Lines are drawn 
"by hand" to make the upper envelope differences clear. Solid 
(dashed) lines trace the upper envelope of the thin-disk (thick- 
disk) star data. The dotted lines trace the upper envelope of thin- 
disk stars older than 2, 5, and 8 Gyr, respectively. Upper limits 
to lithium abundances based on the non detection of the lithium 
6708 A doublet are not shown. 

ure[4]may be influenced (likely dominated) by systematic 
differences in mass, metallicity, and age and, therefore, 
lithium depletions between the two samples. 

Although the thick disk overlaps the thin disk in metal- 
licity, it is systematically older than the thin disk and, 
therefore, the higher mass stars are found exclusively 
in the thin disk. For our first detailed comparison, 
we restrict the thin disk sample to stars with masses 
lower than 1.1 Mq and then plot separately thin-disk 
stars older than 2, 5, and 8 Gyr (Figure [SJ. Previous 
studies have found that most thick-disk stars are older 
than 8 Gyr (|Allende Prieto et al.l 120061 and references 
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Figure 6. Lithium abundance versus age for thin-disk (crosses) 
and thick-disk (circles) stars of mass 0.9 < M/Mq < 1.1. Upper 
limits to lithium abundances based on the non detection of the 
lithium 6708 A doublet are not shown. 

therein). The solid line in each of the four panels of 
Figure [5] denotes the upper envelope to the lithium abun- 
dance for thin-disk stars, i.e., stars younger than 2 Gyr 
fill in the gap between stars plotted in the >2 Gyr panel 
and the solid line. The upper envelope sits progressively 
at lower An for older stars, in fact matching the upper 
envelope for the thick disk for stars older than 8 Gyr. 
Even though this is not entirely clear from an inspection 
of the two bottom panels of Figure [5J we note that the 
few thin-disk stars with An ~ 2.4 are among the most 
massive stars shown in that panel, while the thick-disk 
stars with the highest An have an average mass about 
0.1 Mq lower compared to those thin-disk objects. When 
this subtle mass bias is removed, the transition of the An 
upper envelope from old thin disk to thick disk is very 
smooth (see also below the discussion on Figure |6|) . 

For the three thin disk panels shown in Figure the 
An versus [Fe/H] relations are simil ar. This similarity 
is con sistent with the contention by iSestito fc Randica 
(2005) that lithium depletion is ineffective beyond an 
age of 1-2 Gyr, although this may apply only to the stars 
defining the upper envelope of the An versus [Fe/H] rela- 
tion; Figure |3] suggests that lithium depletes throughout 
the main-sequence life of a 1 Mq, solar-metallicity star. 

Figure [5] shows that over the [Fe/H] range for the thin 
disk, the upper envelope of the An data is almost in- 
dependent of [Fe/H], at a level that decreases only very 
slightly as the age limit is increased from 2 to 5 and then 
8 Gyr. This maximum lithium abundance approaches 
the maximum An shown by thick-disk stars over their 
full [Fe/H] range which extends to lower [Fe/H] than do 
the thin-disk stars. 

A second and related comparison of thin and thick 
disk is presented in Figure |6l which shows lithium abun- 
dance versus age for thin- and thick-disk stars with M = 
1.0 ± 0.1 Mq. The mass interval chosen for this compari- 
son ensures that relatively massive thin-disk stars, which 
have suffered little lithium depletion, and relatively low- 
mass thick-disk stars, which have endured severe lithium 
depletion, are not compared to each other. This mass 
range also contains an important number of both thin- 
and thick-disk stars, allowing a better comparative anal- 
ysis. Figure IH1 confirms the impression from Figure[S]that 
over the age range of thick-disk stars there is very little 
difference in the distribution of lithium abundances for 
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Figure 7. Lithium abundance versus [Fe/H] for our kinematically- 
selected thick-disk (filled circles) and halo stars (open symbols). 
Upper limits to lithium abundances based on the non detection of 
the lithium 6708 A doublet are not shown. 

the two disks. In fact, this thick-disk star An versus 
age trend, particularly at ages older than about 8 Gyr, 
appears to be a natural continuation of the thin disk 
sample. The fact that the few relatively young thick- 
disk stars (age< 8 Gyr) have somewhat low An probably 
owes to the fact that they are among the most massive 
of this group of stars. 

4.2.3. A Primordial Lithium Abundance 

Figure [7] shows the lithium abundances as a function 
of [Fe/H] for our kinematically-selected thick-disk stars 
along with those derived for the halo stars in our sample. 
The maximum An of the most metal-poor stars is com- 
parable with that of the most metal-rich thick-disk stars. 
The star with very high A Li ^ 2.8 at [Fe/H] ~ 0.2 is too 
young to be a "normal" thick-disk member, as argued 
before (see footnote fTTj) . The observation that the maxi- 
mum lithium abundance of thick-disk stars is equivalent 
to that in the more metal-poor stars on the Spite plateau 
recalls the question about whether the Spite plateau rep- 
resents the primordial lithium abundance. If the answer 
is "yes," the puzzle is then to show how the thick disk was 
able to increase its metallicity to [Fe/H] ~ —0.1 with- 
out a significant increase in lithium abundanceQ On 
the other hand, if the answer is "no," the inference that 
the initial lithium abundance was that predicted for the 
WMAP baryon-to-photon ratio (i.e., An — 2.7) leads 
to the conclusion that important surface lithium deple- 
tions have occurred across the entire metallicity and mass 
range spanned by thick-disk and very metal-poor stars. 

In Figure the lithium abundances of thick-disk and 
halo stars from our sample are plotted against their frac- 
tional metallicity Z = Z Q x 10l Fc / H ] (0.638 x 10[ Q / Fc l + 
0.362), where Zq = 0- 14 co rresponds to the Sun, as 
given by lAsplund et al.l (120091). Note that this f ormula 
for Z , which was derived by ISalaris et al.l (|1993D . takes 
into account the impact of a-element enhancement at 
low [Fe/H]. In our case, we assume [a/Fe] = +0.4 for 
[Fe/H] < —0.7 and linearly decreasing [a/Fe] values from 

12 Note, however, that most of the iron in the thick disk was 
produced by Type II supernovae whereas lithium originates from 
spallation in the interstellar medium, which could be a much slower 
process and have no significant production over the relatively short 
period when the thick disk was formed. 
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Figure 8. Lithium abundance as a function of fractional metal- 
licity for the kinematically-selected halo and thick-disk stars of 
our sample. Three mass groups are plotted, as described by the 
lower left legend. Quadratic fits to the lowest (dotted line) and 
highest (dash-dotted line) mass groups are overplotted. along with 
their average (dashed line). The thick and thin solid lines at 
Ajji = 2.72 ± 0.06 correspond to the primordial lithium abundance 
inferred from WMAP observations and standard Big Bang nucle- 
osynthesis. 

+0.4 at [Fe/H] = -0.7 to solar at [Fc/H] = 0. The stars 
in Figure E] have been sorted into three mass groups, as 
described by the legend on the lower left part of the Fig- 
ure. Clearly, in these old stars depends on Z, but 
the relation is mass-dependent. At a given Z, lower mass 
stars have lower A^. 

The mass-dependent Au versus Z relations for thick- 
disk and halo stars can be approximated with quadratic 
fits. Least-squares fits were applied to the lowest and 
highest mass groups, and are shown with the dotted and 
dot-dashed lines in Figure El The dashed line corre- 
sponds to the average of these two quadratic fits, and it 
is a very good approximation to the A^ versus Z relation 
of the intermediate mass group. Interestingly, extrapo- 
lating these relations to Z = leads to a "primordial" 
lithium abundance A^ ~ 2.73, a value that is in excellent 
agreement with that derived from WMAP observations 
and s tandard Big Ban g nucleosynthesis Au = 2.72±0.06 
(e.g.. ICoc etaLl 120121) . 

Thus, the data for old thick-disk and halo stars re- 
veals a lithium-mass-metallicity dependence that hints at 
a solution to the primordial lithium abundance problem. 
We note, however, that the physics behind this observa- 
tion is likely related to the extent of the convective zone 
and its dependence on the stars' fundamental properties, 
whereas the problem with the Spite plateau requires ad- 
ditional depletion mechanisms, hence different physics. 
The fact that our extrapolation to Z = results in a 
high primordial lithium abundance consistent with the 
cosmological value could therefore be fortuitous. 

Indeed, the Spite plateau is defined mostly by stars 
with Z < 0.002, which are scarce in our work. Since 
nearly all of those objects have Au ~ 2.2, the quadratic 
fits shown in Figure [8] may not apply to them. Per- 
haps the real Ali versus Z relation turns down towards 
lower Au for Z < 0.001, leading to a low derived pri- 
mordial lithium abundance However, if the masses of 



13 In fact, in a very recent work, Nisscn & Schuster (2012), who 
perform a lithium-mass-metallicity fit and extrapolation to Z = 



these very low metallicity stars (in our mass scale) turn 
out to be lower than those of the lowest mass group in 
our sample, they may exhibit a Z-dependence similar 
to that of the higher mass groups, implying a high pri- 
mordial lithium abundance. We encourage future works 
on lithium abundances of large samples of very metal- 
poor stars to examine this relation carefully, including 
the more metal-rich halo and thick-disk stars. 

4.3. Lithium and Planets 

The connection between stellar surface lithium abun- 
dance and presence of planets is the subject of in- 
tense debate. A number of authors claim to have 
found evidence for e nhanced lithium dep l etion in planet- 
hosting stars (e.g.. iChen fc Zhaol I2006L llsraelian et al.1 
2009; IGonzalez et all 12010ft . but others have concluded 
that the sensitivity of lithium abundance to other stellar 
parameters such as age and metallicity are responsible 
for the lithium abundance differences between stars with 
and without planets, which are therefore the product of 
sampl e biases and not related to the presence of planets 
/e.g.. iRvanl 120001 i Luck fc Heiter| [200l iBaumann etall 
2010£ fGherzi et al.ll2010D~ 

There are 165 stars in our sample around which plan- 
ets have been detected. Their lithium abundance pat- 
terns are shown in Figure [9] along with those correspond- 
ing to the rest of our sample stars, i.e., to stars around 
which planets have not yet been detected. In this Figure, 
the comparison sample of non planet-hosts includes stars 
that have not been searched for sub-stellar mass compan- 
ions. The heterogeneous nature of our sample regarding 
planet detection prevents us from making a more consis- 
tent planet/non-planet comparison. Lists of stars with 
and without planets can be found in the literature, but 
all one can say about the non planet-hosts is that they do 
not have planets below a certain mass and orbital period 
combination which depends on the instrumentation and 
observing strategy used by each group of investigators. 
Taking this into consideration, we prefer to compare our 
sample of known planet-hosts to the sample of otherwise 
ordinary field stars, some of which have been searched 
for planets. 

Other than the lack of planet-hosts on the low lithium 
abundance side of the lithium desert (which is discussed 
below), a simple inspection of Figure |9Ka shows that 
the planet-host stars follow a very similar overall lithium 
abundance pattern compared to the sample of field stars. 
Planets are difficult to detect around warm, massive 
stars, which explains the lack planet-hosts on the hot 
T e g and high mass ends of FiguresEKa and[9]Ab. On the 
other hand, the planet-metallicity correlation, a very well 
defined trend at least for short-period massive planets 
around the dwarf stars which dominate ou r sample (e.g., 
iSantos et al.ll2004UFischer fc Valentll2005[ ). is behind the 
separation of samples in Figure [9jAc. It is interesting 
to notice that the planet-host sample defines a lithium 
dip similar to that observed in the Hyades (Figure [9]Ab) 
better than the comparison sample. The fact that the 
planet-host sample has limited spread in [Fe/H] with a 

nearly identical to ours, but with a linear relation, find that 
their fit using metal-rich halo stars and thick-disk stars does not 
re produce well the very metal-poor star lithium abundance data 
bv lMelendez et al.l IpOTO . 
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Figure 9. Lithium abundance as a function of effective temperature (column a), mass (column b), and metallicity (column c) of the star. 
Known exoplanet hosts are shown with filled circles. Crosses correspond to stars for which no sub-stellar mass companions have been found 
yet, including stars which have not been searched for planets. Row A shows the entire sample. Row B shows only stars in the effective 
temperature range from 5650 to 5820 K. In row C, only stars with a mass of M = 1.00 ± 0.05 Mq are shown. Row D shows only stars in 
the effective temperature range from 5820 to 6190 K, i.e., in the T e jr range of the lithium desert. 



mean metallicity that is super-solar and the Hyades have 
[Fe/H] ~ +0.1 while the comparison sample covers a very 
wide [Fe/H] range could be the reason for this observa- 
tion. 

In order to better understand lithium depletion in 
planet hosts, it is customary to restrict the sample to 
a narrow range of effective temperature or mass. This 
practice, however, can be risky. To illustrate this, in Fig- 
ure [9)3 we show the lithium abundance patterns of stars 



in the narrow T e g range from 5650 to 5820 K. A simple 
look at Figure 03a suggests that planet-hosts have sys- 
tematically slightly lower lithium abundances. There are 
no planet-hosts with A^i > 2 in this effective tempera- 
ture range, but a good number of comparison stars reside 
there. Closer inspection of this sample, however, allows 
us to better understand this trend and to realize that 
planets are not responsible for it. Figure |9j3b shows that 
the planet-host sample is on average more massive than 
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the comparison sample. Although one could expect these 
stars to suffer less lithium depletion due to their thinner 
convective envelopes, we must notice that these objects, 
on average, are also more metal-rich than the compar- 
ison sample, as shown in Figure |9j3c. The majority of 
non planet-hosting stars with Ai\ > 2 in Figure [9)3 are 
either more massive than M ~ 1.1 Mq or more metal- 
poor than [Fe/H] ~ —0.1. Both conditions imply higher 
lithium abundance, as shown by Figure [2j 

Selecting stars within a narrow mass range does not 
help understanding the connection between lithium and 
planets either. For example, Figure [Spa shows that 
when stars of nearly one-solar mass are selected (M = 
1.00±0.05 Mq), the range of T e g covered by the sample is 
very large, owing to the different ages and metallicities of 
the otherwise "Sun-like" stars. As with the T c g selection, 
this constraint in mass seems to also imply enhanced 
lithium depletion in planet-host stars. Figure |HPb would 
suggest that the planet-host sample is shifted downwards 
in Au by about 0.5 dex relative to the comparison stars. 
However, it is clear that this trend is due to the inclu- 
sion of a good number of relatively cool planet hosts with 
relatively massive convective envelopes, which could be 
slightly evolved 1 Mq stars that have had plenty of time 
to deplete surface lithium. In addition, the planet-host 
star sample is again biased towards higher metallicities, 
which implies more lithium depletion at any given stellar 
age for stars of about one solar mass, as shown in Fig- 
ure |3l Figure [9]Cc shows in fact a very clear correlation 
between lithium abundance and [Fe/H] for stars of one 
solar mass, a correlation that is independent of whether 
the star has a planet or not. 

Thus, the analysis of our large sample indicates that 
previously claimed patterns of enhanced lithium deple- 
tion in planet-host stars are a pro duct of sample biases in 
stellar mass, age, and metallicity. iBaumann et al.l (|2010l ) 
arrived at the exact same conclusion in their study of 
Sun-like stars. A system of particular interest in this 
context is the binary pair 16 Cyg, which consists of two 
solar analog stars, one of which, the secondary, is known 
to host a gas giant planet and is also lithium poor rela- 
tive to the primary, whi ch appears to not have sub-stellar 
mass companions (e.g.. ICochran et all 11997; Kin g et all 
19971). Using our lit hium abundance catalog, however, 
Ramirez et al.l ()2011[ ) have shown that the lithium abun- 
dance difference seen in the 16 Cyg pair can be fully 
explained by the slightly different masses and metallici- 
ties of the stars, leaving little or no room for a possible 
lithium-planet connection (see their Figure 11). 

The only region of stellar parameter space that seems 
to separate planet hosts and comparison stars into high 
lithium and low lithium abundance samples is the re- 
gion of the lithium desert (Figure [9j3a) . Interestingly, 
in this region (5820 < T c g < 6190 K), the planet-hosts 
seem to have higher lithium abundance than ordinary 
field stars. More specifically, there are no planet-hosts 
on the low lithium abundance side of the lithium desert. 
Mass and metallicity effects cannot fully explain what is 
observed here. The mass distributions of the two samples 
are very similar (Figure 03b) and, even though planet- 
hosts are, as usual, more metal-rich than the comparison 
stars, in the region of [Fe/H] where there is significant 
overlap (—0.1 < [Fe/H] < +0.1), one can still see that 
there are comparison stars with low lithium abundances 



(Figure [5pc) . Comparison stars with low lithium abun- 
dances in the desert region are still present if we restrict 
the sample to [Fe/H] > -0.2 and M < 1.2 Mq, i.e., 
when we make the two samples as similar as possible in 
their atmospheric parameters. A number of them (six) 
still remain there if we further restrict the comparison 
sample to stars that have been searched for planets but 
none have been detected around them yet. As explained 
in Sect. 14.11 stars on the low lithium abundance side of 
the desert are older than about 2 Gyr. Only one of the 
planet-hosts in the T e g range of the desert is younger 
than 2 Gyr. Thus, the fact that the planet-host stars in 
Figure [HP have lithium abundances at least about one 
order of magnitude larger than those on the low lithium 
abundance side of the desert appears not to be due to an 
age bias either. 

Our data suggest that some mechanism of planet- 
star interaction (not necessarily related to planet forma- 
tion) prevents the star from undergoing the sudden sur- 
face lithium abundance drop responsible for the lithium 
desert. Nevertheless, we must practice extreme caution 
with this type of statement. Stars on the low side of 
the desert are relatively rare. Also, most of the planets 
discovered around our sample stars are gas giants. The 
existence of smaller planets around the stars on the low 
lithium abundance side of the lithium desert must not 
be ruled out. A systematic search for planets around the 
stars in the lithium desert region would be interesting to 
determine whether or not lithium abundances can really 
discriminate planet-hosts from non planet-host stars, at 
least within this narrow T c g window. 

5. CONCLUSIONS 

We have measured stellar parameters and lithium 
abundances of 671 stars, of which, to the best of our 
knowledge, 319 have their lithium abundances derived 
for the first time. Then, using data from the litera- 
ture, we have constructed a catalog of stellar parameters, 
lithium abundances, masses, and ages for 1381 nearby 
FGK dwarf and subgiant stars. These data are used 
to investigate stellar lithium depletion inside stars, the 
Galactic chemical evolution of lithium, and the impact of 
planet formation on the internal evolution of cool stars. 

A lithium desert is observed in the versus T e ff plane 
from about T c g — 5950 to 6100 K. Stars are found above 
and below Au ~ 1.8, creating a gap of about 0.5 dex in 
A^i. Many of the stars on the low lithium abundance side 
of the desert have only upper limits to A^, implying that 
the gap could be deeper. Detailed inspection of Au ver- 
sus mass, age, and metallicity plots suggests that a short- 
lived process depletes lithium on stellar surfaces during 
their main-sequence or subgiant phases, at least for stars 
with masses between 1.1 and 1.3 Mq. Stars from the 
lithium dip, a feature observed in some young open clus- 
ters, which have now evolved and separated from their 
parent clusters, may contribute to the low A^i side of 
the lithium desert, but they are unlikely to be the main 
source of stars in that region of the versus T e ff plane. 

We have found evidence of lithium depletion on the 
main sequence for stars other than the Sun, and have 
revealed a dependence of lithium depletion on the stel- 
lar metallicity, which is, however, not straightforward. 
Metal-rich stars deplete lithium quicker than metal-poor 
stars for M < 1.1 Mq, but the opposite is true for higher 
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masses. Time-scales of main-sequence lithium depletion 
depend both on mass and metallicity, but not in a simple 
way. 

We have attempted to detect, for the first time to the 
best of our knowledge, a difference in lithium abundances 
for thin- and thick-disk stars. However, the stars with 
maximum lithium abundance in each of these groups, i.e., 
those with the least amount of lithium depletion, have 
different fundamental properties. Thin-disk stars tend 
to be younger and more massive than thick-disk stars, 
leading to thin/thick disk lithium abundance differences 
that reflect different degrees of lithium depiction rather 
than differences in lithium enrichment of the interstel- 
lar medium. Nevertheless, we find that the maximum 
lithium abundance of thick-disk stars is very similar to 
the Spite plateau value defined by very old halo stars, 
suggesting that it may extend up to [Fe/H] = —0.1. 

Combining thick-disk and halo stars, we find a well- 
defined relation between lithium abundance and frac- 
tional metallicity, which is mass-dependent. Extrapo- 
lated to zero metallicity, this relation implies a primor- 
dial lithium abundance of Au = 2.73, i.e., a value that is 
consistent with WMAP observations and standard Big 
Bang nucleosynthesis. Although this hints at a solution 
to the primordial lithium abundance problem, we note 
that our work does not include many very metal-poor 
stars, which could affect the shape of the lithium-mass- 
metallicity relation, and therefore the primordial lithium 
abundance inferred from its extrapolation to zero metal- 
licity. 

Our investigation of lithium depletion in planet-host 
stars relative to ordinary field stars shows that a pat- 
tern of enhanced lithium depletion in planet-hosts can be 
created by restricting samples to narrow ranges of T e g 
or mass. However, age and metallicity effects must be 
taken into account before drawing any conclusions from 
this biased pattern. When all relevant stellar parame- 
ters are properly taken into account, we find no evidence 
for enhanced lithium depletion in planet-hosts. In fact, 
we find instead that in the region of the lithium desert, 
no planet hosts are found with low lithium abundances. 
Only a dedicated search for planets around stars in the 
lower side of the lithium desert can tell us whether this 
is something related to a planet-star interaction or if it 
is due to sample biases, but with the data in hand the 
former scenario is favored. 

I.R.'s work was performed under contract with the 
California Institute of Technology (Caltech) funded by 
NASA through the Sagan Fellowship Program. J.R.F. 
acknowledges support from the summer student program 
at Carnegie Observatories. D.L.L. wishes to thank the 
Robert A. Welch Foundation of Houston, Texas for sup- 
port through grant F-634. The authors thank the ref- 
eree, Dr. Luca Sbordone, for his suggestions to improve 
our manuscript. 

REFERENCES 

Alibes, A., Labay, J., & Canal, R. 2002, ApJ, 571, 326 

Allcndc Pricto, C, Barklcm, P. S., Lambert, D. L., & Cunha, K. 

2004, A&A, 420, 183 
Allende Prieto, C, Beers, T. C, Wilhelm, R., Newberg, H. J., 

Rockosi, C. M., Yanny, B., & Lee, Y. S. 2006, ApJ, 636, 804 



Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, 

ARA&A, 47, 481 
Bagnulo, S., Jchin, E., Ledoux, C, Cabanac, R., Melo, C, 

Gilmozzi, R., & The ESO Paranal Science Operations Team. 

2003, The Messenger, 114, 10 
Balachandran, S. 1990, ApJ, 354, 310 
— . 1995, ApJ, 446, 203 

Balachandran, S. C, Mallik, S. V., & Lambert, D. L. 2011, 

MNRAS, 410, 2526 
Barnes, S. A. 2007, ApJ, 669, 1167 

Baumann, P., Ramirez, I., Melendez, J., Asplund, M., & Lind, K. 

2010, A&A, 519, A87 
Bensby, T., Feltzing, S., Lundstrom, I., & Ilyin, I. 2005, A&A, 

433, 185 

Boesgaard, A. M., Budge, K. G., & Ramsay, M. E. 1988, ApJ, 
327, 389 

Boesgaard, A. M., & Tripicco, M. J. 1986, ApJ, 302, L49 
Bouvicr, J. 2008, A&A, 489, L53 

Brugamyer, E., Dodson-Robinson, S. E., Cochran, W. D., & 

Sneden, C. 2011, ApJ, 738, 97 
Casagrandc, L., Ramirez, I., Melendez, J., Bessell, M., & 

Asplund, M. 2010, A&A, 512, 54 
Casagrande, L., Schonrich, R., Asplund, M., Cassisi, S., Ramirez, 

I., Melendez, J., Bensby, T., & Feltzing, S. 2011, A&A, 530, 

A138 

Castro, M., Vauclair, S., Richard, O., & Santos, N. C. 2008, 

Mem. Soc. Astron. Italiana, 79, 679 
Cayrel de Strobel, G., & Cayrel, R. 1989, A&A, 218, L9 
Chaname, J., & Ramirez, I. 2012, ApJ, 746, 102 
Charbonnel, C, & Talon, S. 2005, Science, 309, 2189 
Chen, Y. Q., Nissen, P. E., Benoni, T., & Zhao, G. 2001, A&A, 

371, 943 

Chen, Y. Q., & Zhao, G. 2006, AJ, 131, 1816 
Coc, A., Goricly, S., Xu, Y., Saimpert, M., & Vangioni, E. 2012, 
ApJ, 744, 158 

Cochran, W. D., Hatzes, A. P., Butler, R. P., & Marcy, G. W. 

1997, ApJ, 483, 457 
Cutri, R. M., et al. 2003, 2MASS All Sky Catalog of point 

sources. (The IRSA 2MASS All-Sky Point Source Catalog, 

NASA/IPAC Infrared Science 

Archive. http://irsa.ipac.caltech.edu/applications/Gator/) 
D'Antona, F., & Mazzitelli, I. 1994, ApJS, 90, 467 
Fields, B. D. 2011, Annual Review of Nuclear and Particle 

Science, 61, 47 
Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102 
Gaidos, E. J., Henry, G. W., & Henry, S. M. 2000, AJ, 120, 1006 
Ghezzi, L., Cunha, K., Smith, V. V., & de la Reza, R. 2010, ApJ, 

724, 154 

Gilmore, G., & Reid, N. 1983, MNRAS, 202, 1025 
Gonzalez, G. 2008, MNRAS, 386, 928 

Gonzalez, G., Carlson, M. K., & Tobin, R. W. 2010, MNRAS, 
403, 1368 

Isaacson, H., & Fischer, D. 2010, ApJ, 725, 875 

Israelian, G., et al. 2009, Nature, 462, 189 

Kim, Y., Demarque, P., Yi, S. K., & Alexander, D. R. 2002, 

ApJS, 143, 499 
King, J. R., Dcliyannis, C. P., Hiltgcn, D. D., Stephens, A., 

Cunha, K., & Boesgaard, A. M. 1997, AJ, 113, 1871 
Lambert, D. L., & Reddy, B. E. 2004, MNRAS, 349, 757 
Lind, K., Asplund, M., & Barklem, P. S. 2009, A&A, 503, 541 
Luck, R. E., & Heiter, U. 2006, AJ, 131, 3069 
Melendez, J., Bergemann, M., Cohen, J. G., & et al. 2012, A&A, 

in press 

Melendez, J., Casagrande, L., Ramirez, I., Asplund, M., & 

Schuster, W. J. 2010, A&A, 515, L3 
Mendes, L. T. S., D'Antona, F., & Mazzitelli, I. 1999, A&A, 341, 

174 

Mermilliod, J., Mermilliod, M., & Hauck, B. 1997, A&AS, 124, 
349 

Michaud, G. 1986, ApJ, 302, 650 

Mishcnina, T. V., Soubiran, C, Kovtyukh, V. V., & Korotin, 

S. A. 2004, A&A, 418, 551 
Nissen, P. E., & Schuster, W. J. 2012, A&A, in press 

(arXiv/1205.4122) 
Nordstrom, B., et al. 2004, A&A, 418, 989 
Perryman, M. A. C, et al. 1997, A&A, 323, L49 



LITHIUM ABUNDANCES IN NEARBY FGK DWARF AND SUBGIANT STARS 



15 



Pinsonneault, M. H. 2010, in IAU Symposium, Vol. 268, IAU 
Symposium, cd. C. Charbonnel, M. Tosi, F. Primas, & 
C. Chiappini, 375-380 

Pinsonneault, M. H., DePoy, D. L., & Coffee, M. 2001, ApJ, 556, 
L59 

Pizzolato, N., Maggio, A., Micela, G., Sciortino, S., & Ventura, P. 

2003, A&A, 397, 147 
Ramirez, I., Allcnde Prieto, C, & Lambert, D. L. 2007, A&A, 

465, 271 

Ramirez, I., Melendez, J., Cornejo, D., Roederer, I. U., & Fish, 

J. R. 2011, ApJ, 740, 76 
Reddy, B. E., Lambert, D. L., & Allende Prieto, C. 2006, 

MNRAS, 367, 1329 
Reddy, B. E., Tomkin, J., Lambert, D. L., & Allende Prieto, C. 

2003, MNRAS, 340, 304 
Richard, O., Michaud, G., & Richer, J. 2005, ApJ, 619, 538 
Romano, D., Matteucci, F., Molaro, P., & Bonifacio, P. 1999, 

A&A, 352, 117 
Ryan, S. G. 2000, MNRAS, 316, L35 

Ryan, S. G., Kajino, T., Beers, T. C, Suzuki, T. K., Romano, D., 

Matteucci, F., & Rosolankova, K. 2001, ApJ, 549, 55 
Salaris, M., Chieffi, A., & Straniero, O. 1993, ApJ, 414, 580 
Santos, N. C, Israelian, G, & Mayor, M. 2004, A&A, 415, 1153 
Scholz, A., Coffey, J., Brandekcr, A., & Jayawardhana, R. 2007, 

ApJ, 662, 1254 
Sestito, P., & Randich, S. 2005, A&A, 442, 615 



Sncdcn, C. A. 1973, PhD thesis, The University of Texas at 
Austin 

Soderblom, D. R., Fedele, S. B., Jones, B. F., Stauffer, J. R., & 

Prosscr, C. F. 1993, AJ, 106, 1080 
Soubiran, C. 1993, A&A, 274, 181 
Sousa, S. G., et al. 2008, A&A, 487, 373 
Spite, F., & Spite, M. 1982, A&A, 115, 357 
Steigman, G. 2010, in IAU Symposium, Vol. 265, IAU 

Symposium, ed. K. Cunha, M. Spite, & B. Barbuy, 15-22 
Strassmcicr, K., Washucttl, A., Granzer, T., Scheck, M., & 

Weber, M. 2000, A&AS, 142, 275 
Takcda, Y., Honda, S., Kawanomoto, S., Ando, H., & Sakurai, T. 

2010, A&A, 515, A93 
Talon, S., & Charbonnel, C. 1998, A&A, 335, 959 
Theado, S., & Vauclair, S. 2012, ApJ, 744, 123 
Travaglio, C, Randich, S., Galli, D., Lattanzio, J., Elliott, L. M., 

Forestini, M., & Ferrini, F. 2001, ApJ, 559, 909 
Tull, R. G., MacQueen, P. J., Sneden, C, & Lambert, D. L. 1995, 

PASP, 107, 251 
Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 
van Leeuwen, F. 2007, A&A, 474, 653 
Wright, J. T., et al. 2011, PASP, 123, 412 

Yi, S., Demarque, P., Kim, Y., Lee, Y., Ree, C. H., Lejcune, T., 
& Barnes, S. 2001, ApJS, 136, 417 



